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Some characteristics of midvelocity emissions in semiperipheral heavy-ion collisions at Fermi ener- 
gies are discussed in the framework of a multifragmenting scenario. We report on binary dissipative 
collisions of "^Nb + ^^Nb at 38^ MeV in which we measured an abundant emission of particles and 
fragments not originated from the usual evaporative decay of hot primary fragments. We checked 
the compatibility of these emissions with the multifragmentation of a source which forms in the 
overlap region. One can fairly well reproduce the data assuming a hot and dilute source, possibly 
more n-rich than the initial nuclei; the results appear to be insensitive to the source size. 

PACS numbers: 25.70.Lm, 25.70.Pq 
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I. INTRODUCTION 

The recent literature on experimental heavy-ions 
physics in (semi-) peripheral collisions at Fermi energies 
(see e.g. [H, 0, H, |j, [EIIjH, Hi) often focuses on the emis- 
sion of particles and fragments from the phase-space re- 
gion in between the reacting nuclei, called midvelocity 
emission. This phenomenon has been already investi- 
gated by several authors, as a function of the size of the 
system and of the impact parameter, but its origin is far 
from being understood; in particular, it is not clear if one 
deals with only one hot source, if the source(s) is(are) 
locally thermalized and what are the regions accessed 
by these systems in the equation-of-state phase diagram. 
From a theoretical point of view, understanding the mid- 
velocity emissions is linked to the more general issue of 
isospin dynamics |T3| , which are supposed to be ruled 
by the symmetry-energy term. This term is unknown 
at the low barion densities predicted for the midveloc- 
ity matter when a neck develops between the interacting 
nuclei during the fast separation phase 0, [U . 

Several experiments (e.g. [1, [l2|) have shown that 
the characteristics of the midvelocity emissions are very 
different with respect to the emissions from the quasi- 
projectile {QP) [and quasi-target (QT)] and standard 
evaporation codes can give a good reproduction of the 
QP emissions, but not of the midvelocity emissions. 
From an accurate energy balance [3|] it was found that 
the QP and QT become more and more excited with de- 
creasing impact parameter. At the same time, already 
for peripheral impacts, the energy deposited in the mat- 
ter emitted at midvelocity accounts for a large fraction of 
the energy globally dissipated into the internal degrees of 
freedom; it is possible that the resulting large excitations 
easily overcome the limit of 3-4^ MeV for which many 
authors have shown the onset of instabilities in nuclei. 



Therefore, we found appropriate to verify the compati- 
bility of the experimental results with calculations per- 
formed in the frame of the Statistical Multifragmentation 
Model (SMM) [I3|,llj,|il. 



II. CHARACTERIZATION OF THE REACTION 
SCENARIO 



The experimental data used for the comparison with 
the SMM calculations are those of the ^^Nb -\- ^^Nb col- 
lision at 38^ MeV [6]. The experimental setup, the anal- 
ysis procedures and the method for separating the mid- 
velocity emissions from the evaporative decays are de- 
scribed in detail in Refs. [1, U^. Here we briefly remind 
that binary events were fully characterized by detecting 
in coincidence the QP and the QT, together with the as- 
sociated emissions of light charged particles (LCP, Z < 2) 
and intermediate mass fragments (IMF, 3 < Z < 7); iso- 
topic discrimination was obtained for Z=l and charge 
identification for reaction products up to Z—7. Differ- 
ent aspects of the midvelocity emissions were presented 
elsewhere [H, S [U [1-8] and their characteristics were com- 
pared with those of the sequential decay of the two main 
reaction products, the excited QP and QT. 

All the observables related to the hot QP and QT are 
compatible with the statistical decay of a source at nor- 
mal density, which can reach excitation energies up to 
about 3^ MeV (T < 6 MeV) Wi^L in agreement with 
the findings of other authors [4. Il2l. [l9l . l20l [2l| when 
comparable windows of impact parameter are selected. 
In particular, the LCP multiplicities and their ratios, as 
well as the neutron content of hydrogen isotopes (see for 
instance Figs. 10-12 of Ref. ^6]), are in good agree- 
ment with the predictions of commonly used evaporation 
codes. (At variance with other authors, e.g. Ref. [12], 
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FIG. 1: (color online) Experimental observables for mid- 
velocity (solid circles) and QP (open circles) emissions as a 
function of TKEL: (a) Neutron to proton ratio for hydro- 
gen isotopes; (b) Ratio of the yield of protons to a-particles. 
The lines through the experimental points are guides to the 
eye. The vertical bars on the right of the experimental points 
correspond to the range of acceptable solutions of SMM cal- 
culations for the midvelocity source (solid bars) and the QP 
(open bars), as explained in the text. 



we don't need to invoke a statistical multifragmentation 
mechanism for the observed decay of QP and QT). 

On the contrary, the midvelocity emissions present fea- 
tures which make them strongly different from the usual 
decay of hot rotating nuclei. This is clearly shown by 
the results of Fig. [TJ Part (a) of the figure presents the 
average experimental isospin-ratio of hydrogen, {N/Z)jj, 
both for the midvelocity emissions (solid circles) and for 
the QP decay (open circles), as taken from Fig. 12 of 
Ref. . Part (b) of Fig. [T] presents the experimental 
ratio of the proton/a-particle yields, (Yp/Ya)^ again for 
the midvelocity emissions (solid circles) and for the QP 
decay (open circles) ; this part has been derived from the 
multiplicities of Fig. 5 of the same Ref. @. 

The experimental results are presented as a function of 
a variable, TKEL, which is used here just as an "order- 
ing parameter" to classify the events in bins of decreasing 
impact parameter ^ . This use of the variable TKEL has 



^ TKEL is defined as the difference between the initial center-of- 
mass energy EjJ; ™' of the collision and the cm. kinetic energy 
of the outgoing QP and QT, assuming binary kinematics 



been extensively discussed in a previous paper Q , where 
it was shown that - even at Fermi energies - TKEL is 
strictly correlated with the relative velocity between QP 
and QT and, as such, it is an indicator of the violence of 
the collision (and hence of the impact parameter); this 
interpretation has been recently adopted also in a theo- 
retical description of the reaction dynamics [ll| . The es- 
timated correspondence of TKEL with impact parameter 
for the data of this paper is displayed by the arrow at the 
top margin of Fig. [TJ It shows that the presented data re- 
fer to semi-peripheral collisions in the impact-parameter 
range between 8 and 11 fm. More central events have 
not been used because, with decreasing relative velocity, 
the separation between midvelocity emissions and QP 
(or QT) evaporation becomes more and more uncertain 
and hence the interpretation of the results would become 
less and less reliable. 

From the data of Fig. [T] it is apparent that the mid- 
velocity emissions present some substantially different 
characteristics with respect to the usual evaporative pro- 
cesses. While the {N/Z)h ratio for particles emitted by 
the QP (open circles) is rather low and it was found 
compatible with the results of Gemini calculations for 
the decay of a ^^Nb nucleus at normal density with the 
appropriate excitation energy [6(], the same ratio for the 
midvelocity emissions (solid circles) is much higher and 
cannot be reproduced by usual statistical model calcu- 
lations. Even more striking is the fact that the relative 
abundances of the emitted particles are different, display- 
ing some remarkable inversions ; in particular Fig. [Ifb) 
shows that the relative abundances of protons and a par- 
ticles are reversed when going from the evaporative to the 
midvelocity emissions. 

The exact nature of the midvelocity emissions is yet 
not fully ascertained. On the one hand, one observes a 
component which is evidently Coulomb-related [l3,[13| to 
the main reaction products: this suggests that they are 
driven by the reaction dynamics towards rather elongated 
shapes, which quickly decay before complete decoupling 
from the formation stage |l|, [lO]. On the other hand, 
there are the very 'central' emissions, located at small 
velocities in the cm. system: they may be attributed 
to the disassembly of the overlap region (possibly char- 
acterized by high excitations and small densities), or to 
the snap-off of distended neck-like structures formed just 
in between the two separating main products fl] . In any 
case, our previous results [1, [l3| suggest that the energy 
density in this region may reach values well inside the 
region in which nuclear multifragmentation signals have 
been detected. Dynamical codes predict the formation of 
elongated structures in fast dissipative collisions in the 
Fermi energy domain ^, [H [13, El . As to the timescale 
of the formation/emission processes, dynamical calcula- 
tions [§, [ni, 111] for semiperipheral reactions suggest that 
it should be short (about 100-200 fm/c), but still long 
enough to allow the system to reach at least a partial 
equilibration in some degrees-of-freedom at a freeze-out 
stage. The final conclusion about statistical equilibrium 
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can be obtained only with a comprehensive comparison 
with experiments. 

Thus, an hypothetical "source" of the midvelocity 
emissions - supposing that it exists - may have a more 
complex configuration than that considered in the frame 
of a statistical multifragmentation model. Therefore, at 
variance with what one expects in central collisions, in 
peripheral collisions we are aware that SMM calculations 
may not reproduce all details of the midvelocity emis- 
sions, especially concerning the kinematical character- 
istics of fragments, where the role of dynamics is still 
quite important. Nevertheless, the chemical equilibrium, 
responsible for the isospin composition of the fragments 
produced in the source, may be estabhshed. If the model 
is able to capture ~ at least partly ~ some important 
characteristics of the process (which might be, e.g., the 
tendency of some portion of nuclear matter to rapidly 
decay after having reached a partial thermal and chem- 
ical equilibrium, or possibly some kind of dilution be- 
ing subjected to mechanical strains), one may expect 
that the model displays reduced discrepancies with the 
experimental data than the usual low-energy statistical 
models. With this in mind, is seems appropriate to test 
the compatibility of statistical multifragmentation calcu- 
lations with the measured properties of the midvelocity 
emissions. 



III. THE PARAMETERS OF THE SMM 
CALCULATIONS 

The Statistical Muhifragmentation Model (SMM) is a 
well developed code which describes the nuclear disas- 
sembly in various regimes [H, [ij, . Present SMM cal- 
culations include also the decay of the hot fragments pro- 
duced at the freeze-out, via evaporation or Fermi break- 
up [Tsj . Therefore we can directly compare the "final 
state" quantities computed by the code with the experi- 
mental ones, which are by definition secondary, i.e. after 
sequential particle decay. For the comparison with the 
data, we consider here the two experimental variables, 
{N/Z)h and (I^/Yq,), already presented in the previous 
Section: being relative quantities, they are less sensitive 
to systematic uncertainties, both from the experimen- 
tal and theoretical point of view. Other authors used 
other experimental variables for a comparison with SMM. 

In peripheral collisions the input values for the SMM 
calculations are by far less obvious than in the usual ap- 
plications of the model, which mostly refer to a single 
source produced in central collisions [2ll, [2^ . There are 
no reliable hints concerning the appropriate values for 
the source size Ag^ its charge and its excitation en- 
ergy per nucleon e^, as well as their dependence on im- 
pact parameter. In fact, besides the uncertainties due to 
the usual lack of information about the emitted free neu- 
trons, in peripheral collisions the process may be more 
complicated than the simple emission from a source: for 
example, part of the ejected matter could be reabsorbed 



by the QP and QT which, while flying apart, are still 
rather close to each other during the midvelocity emis- 
sion phase ^17|. 

Given these uncertainties on the source parameters, we 
decided to run many SMM calculations with a wide grid 
of input values for Ag^ Zg and eg. A lower limit for the 
charge of the source, Zg^ was given by the experimental 
total charge of the midvelocity emissions, {Zmidv), which 
smoothly increases with TKEL, up to about 12-15 for the 
highest TKEL values considered in this paper. In other 
words, the lower limit for Zg corresponds to the extreme 
assumption that the source completely disassembles into 
the measured total charge {Zmidv)- As an upper limit, 
one might have taken the overlap matter in a participant- 
spectator picture (with or without taking into account 
the diffuseness); however, in order not to bias the results 
with a too restricted choice, we preferred to use in all 
cases a value of 50 charge units, which is rather large (it 
corresponds to more than half of the total charge of our 
system and we recall that we are dealing with periph- 
eral collisions), checking afterwards the sensitivity of the 
results on this parameter. 

The second input parameter is the mass of the source, 
Ag. Multifragmentation processes are much more sen- 
sitive to the neutron content of the source, {N/Z)g = 
{Ag — Zs)/Zg, than to the source size itself. In fact the 
isospin dynamics are supposed to play a relevant role in 
the fragment production, either via the density depen- 
dence of the symmetry energy 0, [2^ or via polarization 
effects of the nuclear matter due to the Coulomb field 
[l^ll^. Our symmetric system has an initial N/Z value 
of 1.27, but we have no other clue on the isospin evolu- 
tion during the interaction and hence on the N/Z values 
appropriate for the various portions of the system. If no 
free midvelocity neutrons were emitted at all, from the 
experimental data one would obtain for the midvelocity 
emissions a low (maybe unrealistic) value of N/Z=1.02- 
1.06 at TKEL >400MeV. Therefore, again in order not 
to bias the results, for each Zg we allowed a large vari- 
ation of the source size Ag, so that the resulting ratio 
{N/Z)g was comprised in the wide range between 0.9 
and 2.0. 

Concerning the excitation of the source, Cs, the ex- 
perimental value obtained via calorimetry - assuming a 
complete disassembly of the midvelocity source into mid- 
velocity emissions, plus some neutrons (Fig. 3 in Ref. Q 
or Fig. 9 in Ref. Q) - points to high excitations for the 
overlapping matter, well above the values deduced for the 
QP and QT . In the calculations we let span the range 
from 4 up to 9 MeV, which approximately corresponds 
to the vaporization limit. We neglect any radial-energy 
contribution to the excitation, as it has been shown that 
the compression-expansion cycle is weak in peripheral re- 
actions Moreover SMM results on fragment yields 
are insensitive to radial flow . 

Another parameter of SMM calculations is the freeze- 
out density of the source. Values of p ranging from po/6 
[1] to po/2> 0, [13, [131 have been commonly used in 
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SMM calculations (with po normal nuclear density). The 
possibility that dilute systems are formed also in semipe- 
ripheral collisions has been suggested by BNV calcula- 
tions ^ which show that, for a system of similar size 
(i24Sn at 35A MeV (b=6 fm)), very low densities 
can be attained in the neck region within 160-250 fm/c, 
especially when a soft equation-of-state is assumed. We 
run SMM with the value p — po/G and we verified that 
the results do not appreciably change when increasing 
the density up to p = po/S. 

A last word on the source shape: the inclusion of de- 
formed shapes for the source modelization is in princi- 
ple relevant. As already said, elongated transient nu- 
clear systems are likely to be formed in semiperiph- 
eral collisions; recently, exotic shapes have been sug- 
gested even for central event sources, thus leading to 
multifragment calculations for non-spherical configura- 
tions '28]. However we think that our data concerning 
isotopic yields of lightest fragments are insufficient to in- 
vestigate with SMM also this degree-of-freedom, which 
therefore has been neglected here (only spherical geome- 
try is assumed). 



IV. COMPARISON WITH THE DATA AND 
DISCUSSION 

Among the many triples {As,Zs,es) defining the 
sources on the grid for the SMM calculations, we retained 
only those for which both calculated values of {N/Z)h 
and (Yp/Ya) differed from the corresponding experimen- 
tal values by less than 25%. Out of over 5000 triples 
used for the input grid, only about 200 were found to 
pass the above mentioned criterion of goodness at each 
TKEL. The values of {N/Z)h and {Yp/Ya) for the so de- 
fined "good" SMM sources are shown, for each TKEL, 
by squares in Fig. [TJ Since all the "good" sources tend 
to give values of these two parameters quite close to each 
other, the squares tend to bunch up into the vertical solid 
bars which are visible in Fig. [U slightly shifted -for the 
sake of clarity- to the right of each experimental point. 

It is worth noting that many slightly different sources, 
all with rather reasonable parameters, are able to simul- 
taneously reproduce both experimental correlations of 
Fig. [Tl including their dependence on impact parame- 
ter; this is a result which was not obvious a priori. 

As it was already shown in Ref. the midvelocity 
emissions are appreciably richer in deuterons and tritons 
(with respect to the protons) than the evaporative emis- 
sions, a feature which is reproduced by the SMM cal- 
culations. The results of the SMM calculations tend to 
be even more n-rich than the experimental data, but the 
discrepancy is not very large ^ . Also the inversion in the 



relative abundances of protons and a-particles between 
midvelocity and evaporative data, which is displayed by 
the ratio (Xp/Ya) in Fig. [TJa), is well reproduced by the 
same "good" SMM calculations. 

This comparison shows that the examined features of 
midvelocity emissions bear a closer resemblance with the 
expectations of a multifragmentation model than with 
those of a sequential decay. To give more support to this 
last statement, we note that, in its recent versions, the 
SMM code has been upgraded to treat also the 'stan- 
dard' evaporation of hot nuclei at moderate excitations 
and normal density. Therefore we have run the SMM 
code also in this mode (with input parameters similar to 
those used for previous GEMINI calculations, see Fig. 12 
of Ref. 0) and also this comparison with the experimen- 
tal data for the QP decay is presented by the open bars 
to the right of the open circles in Fig. [BJa). It is worth 
noting that the large difference in {N /Z)h between evap- 
orative and midvelocity emissions is well reproduced by 
SMM. We did not try to calculate {Yp/Ya) for the evap- 
oration of the QP within SMM, since this ratio strongly 
depends on the angular momentum of the emitting source 
[30, |3l| . This would require the determination of angular 
momentum as a function of the impact parameter for the 
QP source and it would divert from the main subject of 
the paper, which concerns the midvelocity source. 

It may be worth examining in more detail the param- 
eters of the "good" sources which lead to the observed 
agreement with the experimental points for the midve- 
locity emissions (circles in Fig. [1]). For the source size 
As, rather flat distributions are obtained, thus indicat- 
ing the insensitivity of the presently considered isospin 
quantities to this SMM parameter; a similar insensitivity 
of other experimental quantities to the source size was 
already noted in Refs. 0, 3|. Concerning the energy 
density, es, in all cases the best agreement is obtained 
for rather large excitation energies, well inside the mul- 
tifragmentation region: most good triples have eg in the 
range 7~9 MeV, with a very weak tendency to increase 
with increasing TKEL. 

Coming to the neutron content, the {N/Z)s of the 
"good" sources is presented as a function of TKEL in 
Fig. [21 where all distributions have been normalized to 
unity. It is interesting to note that the average value of 
{N/Z)s (indicated by the arrows in Fig. ^ displays a 
definite increasing trend with increasing TKEL. Only in 
the most peripheral bin it remains below that of the col- 
liding system (1.268, indicated in figure by the vertical 
line), while for the other less peripheral bins it rises and 
reaches values around and above 1.4. If the source size 
increases with decreasing impact parameter, this might 
naturally tend to privilege higher N/Z values. To ex- 



^ We note that the midvelocity yields, obtained in Q by subtract- 
ing the evaporative component from the total, actually refer to all 



not-equilibrated emissions, including those which are anisotrop- 
ically distributed on the Coulomb ridges of the QP and QT 
[TtI. I29II and might be less exotic in their neutron content. 
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FIG. 2: {N/Z)s distributions for the midvelocity multifrag- 
menting sources of the SMM calculations which produce the 
results shown in Fig. [1] From top to bottom, the TKEL win- 
dows go from very peripheral to mid-peripheral events. Each 
distribution is normalized to unit area and the arrow shows 
its average value; the vertical line is the N/Z of the system. 



elude this trivial explanation, we checked that no ap- 
preciable correlation exists between {N/Z)s and As; in 
addition, limiting the accepted source sizes to 40 amu 
does not change the trend observed in Fig. [2l One could 
have attributed this trend to the required agreement with 
{N/Z)h, an experimental observable which strongly de- 
pends on deuterons: because of their more diffuse wave 
function, deuteron production could be enhanced in too 
diluted configurations 0. However, it was verified that 
increasing the freeze-out density to p = Po/S produces 
negligible effects on the behaviour of {N/Z)s- 

The question of a possible neutron enrichment is widely 
discussed in the literature. Theoretically, Landau- Vlasov 
calculations with different chemical potentials for neu- 
trons and protons predict -for the neck matter- a neutron 
enrichment with respect to the initial value when an asy- 
stiff equation-of-state is used [1, [Ulllll. For example, in 
a recent work several effects contributing to isospin 
dynamics in semiperipheral reactions were studied. The 
predicted n-enrichment of the neck-matter occurs, even 
for symmetric systems, via the socalled "isospin migra- 
tion" process, which sets in because of the density gra- 
dient between the QP (and QT) and the more diluted 
neck matter. The observed rising trend of [N / Z)s with 



increasing centrality is in agreement with other exper- 
imental indications [T^ . The rather peripheral impact 
parameters addressed in this work (6 « 8-11 fm) do not 
allow to check whether the N/Z of the midvelocity mat- 
ter decreases for substantially more central collisions, as 
it is predicted by some dynamical calculations (see, e.g., 

Experimental results are not easy to compare with each 
other and they usually give indications based on the iso- 
topic analysis of only some species (IMFs or LCPs) emit- 
ted by the multifragmentation of the source. These par- 
tial results are usually interpreted as indications that, 
already in peripheral collisions, the midvelocity source is 
neutron enriched with respect to the the QP source (see 
e-g- 0) 0,11, 113 )• Complete measurements of all kinds of 
products emitted at midvelocity are rare. Combining the 
data of two experiments, one about the free neutrons and 
the other about the charged products, it was concluded 
that globally all the material found in the midveloc- 
ity region is likely to have the same N/Z ratio as the bulk 
matter. However, more recently, a simultaneous isotopic 
analysis of all emitted products (IMFs, LCPs and neu- 
trons) has been performed [l^, leading to the opposite 
conclusion: although the errors are very large, it gives 
support to the idea that the midvelocity matter is more 
n-rich than the initial system and it tends to become even 
more n-rich with decreasing impact parameter. 

In the present case, the results of Fig. [2] are just a 
model-dependent indication that some of the features 
observed at midvelocity point to a possible neutron en- 
richment of the neck region, once a statistical multi- 
fragmentation mechanism is assumed to be responsible 
for this kind of emissions. Whether this possible neu- 
tron enrichment is mainly produced by density depen- 
dent effect]^, [13, or by a lowering of the symmetry 
energy in the dilute hot multifragmenting systems [3^, 
or by some other exotic mechanism, is still an open and 
widely debated question, which deserves further investi- 
gation. 



V. CONCLUSIONS 

In semiperipheral collisions (6 >8 fm) at Fermi ener- 
gies it was found that the experimentally observed in- 
tense emission of reaction products from the midveloc- 
ity region displays some characteristics which are quite 
different from those of the usual evaporation from hot 
nuclei. The excitation energy per nucleon in the over- 
lap region has been estimated to be well above the com- 
monly accepted limit for multifragmentation. Therefore, 
we have investigated in how far the behavior of the above 
mentioned experimental quantities may be reproduced by 
calculations in the frame of the Statistical Multifragmen- 
tation Model (SMM). The main outcome of this investi- 
gation is that it is indeed possible to find source param- 
eters such that the final multifragmentation products re- 
produce reasonably well the observed peculiar features of 
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midvelocity emissions. In particular they reproduce the 
neutron enrichment of the emitted hydrogen isotopes, as 
well as the inversion of relative abundances of protons 
and a-particles with respect to an evaporative process. 
These SMM sources, while showing no preference for a 
particular size Ag, are characterized by rather high exci- 
tation energies (7-8 MeV per nucleon) and by a tendency 
to become more neutron-rich than the colliding system 



with decreasing impact parameter. 
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